arXiv:1504.05246v2 [astro-ph.SR] 12 Apr 2016 


Mon. Not. R. Astron. Soc. 000.[T]-?? (2016) 


Printed 13 April 2016 (MN MTeX style file v2.2) 


Cyclic variations in O—C diagrams of field RR Lyrae stars 
as a result of LITE 


J. Liska^*, M. Skarka^’^, M. Zejda^ Z. Mikulasek^, & S. N. de Villiers^ 

^Department of Theoretical Physics and Astrophysics, Masaryk University, Kotldfskd 2, 61137 Brno, Czech Republic 

^Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences, H-1121 Budapest, Konkoly Thege Mikids 
ut 15-17, Hungary 

^Private Observatory, 61 Dick Burton Road, Plumstead, 7800 Cape Town, South Africa 


Received ...,2016; 


ABSTRACT 

This paper presents an extensive overview of known and proposed RR Lyrae stars in binaries. 
The aim is to revise and extend the list with new Galactic field systems. We utilized maxima 
timings for eleven RRab type stars with suspicious behaviour from the GEOS database, and 
determined maxima timings from data of sky surveys and our own observations. This signif¬ 
icantly extended the number of suitable maxima timings. We modelled the proposed Light 
Time Effect (LiTE) in 0—C diagrams to determine orbital parameters for these systems. In 
contrast to recent studies, our analysis focused on decades-long periods instead of periods in 
the order of years. Secondary components were found to be predominantly low-mass objects. 
However, for RZ Get and AT Ser the mass of the suspected companion of more than one so¬ 
lar mass suggests that it is a massive white dwarf, a neutron star or even a black hole. We 
found that the semi-major axes of the proposed orbits are between 1 and 20 au. Because the 
studied stars belong to the closest RR Lyraes, maximal angular distances between compo¬ 
nents during orbit should at least be between 1 and 13 mas and this improves the chance to 
detect both stars using current telescopes. However, our interpretation of the 0—C diagrams 
as a consequence of the LiTE should be considered as preliminary without reliable spectro¬ 
scopic measurements. On the other hand our models give a prediction of the period and radial 
velocity evolution which should be sufficient for plausible proof of binarity. 

Key words: methods: data analysis - techniques: photometric - stars: horizontal branch - 
stars: variables: RR Lyrae - stars: binaries: general 


1 INTRODUCTION 

RR Lyrae stars play a crucial role in many astrophysical disci¬ 
plines due to their unique characteristics. Since they are evolved, 
luminous, population II stars with an easily recognizable light 
curve, they are used fo r the testing of evolutio nary models (e.g. 
iLee & Demarguel Il990l : IPietrinferni et al.l l2006h and research on 
the evo lution, dynamics, and formation history of the Milky Way 
galaxy ( ICatelaiill2009l : IPrake et al.ll2013l : IStetson et ^12014) . Be¬ 
cause they are detectable over long distances, in combination with 
a known relationship between metallicity, period, and luminos¬ 
ity the y represent suitable distance indicators also at e xtragalactic 
scales dCatelan & Cortes! [20081 : iFiorentino et'^l2010l) . However, 
the location on the horizontal branch (HB) and the luminosity of 
a particular RR Lyrae star also depend on its mass. This parame¬ 
ter has been estimated only indirectly from evolution models and 
pulsation characteristics. Therefore independent determination of 
the mass of these stars is a key task that is sorely needed. If an 
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RR Lyrae would be in a binary system with known orbital param¬ 
eters, its mass could be determined very accurately. Thus finding 
such systems is crucial. 

Unfortunately, only very few RR Lyrae binary stars have been 
unambiguously identified so far (e .g. the field star TU UMa in 
IWade et ^ll999l : rLiska et al.ll2015h . Moreover, the overwhelming 
majority of these stars are only candidates without definite confir¬ 
mation of their binarity. This low number of binaries among RR 
Lyraes is surprising, because RR Lyrae variables are the most nu¬ 
merous class of variables (VSX. Iwatson et ^l2006l) . and recent 
compilations and catalogues with pulsating stars in binaries list 
many hundreds of systems with o th er types of p ulsating stars (e.g. 
ISzatmar^l 199d : ISzabados|[2003al f9: |ZhoiJl20l^ . Well, what could 
be the reason? 

Some possible reasons were discussed, for example, in the 
classical review of RR Lyrae stars provided bv lSmith| ( ll995l) . Since 
RR Lyraes are evolved stars, their possible companion would be in 
the form of a degenerated remnant, another HB star, a red dwarf 
(in the case of main-sequence star), or a substellar object. Thus, 
in most cases, the companion would be significantly less lumi- 
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nous and smaller than the RR Lyrae component itself. They proba¬ 
bly would be only single-line spectroscopic binaries dominated by 
the RR Lyrae component, and possible eclipses w ould be shallow 
and de tectable with difficulty. As it was noted bv ISoszvnski et all 
( l2009h the mass transfer in close binaries would certainly affect 
an RR Lyrae progenitor. Therefore, if we want to discover a binary 
system containing a classical RR Lyrae, i.e. the one that has evolved 
without any interference from a companion, we should search for 
detached systems with wide orbits and with periods longer than 
a few hundreds of days, rather than for short-period binaries. All 
these facts would cause problematic detection due to tiny changes 
over years, or decades. 

Recent research revealed that the situation is not as hope¬ 
less as previous studies suggested. Except for a s trange object re¬ 
sembl ing RR Lyrae-like pulsations, identified byJSosz^nski_et^ 
(201T) in the Ga l actic bulge and studied by IPietrzvnski et al. 
l 2012l) . lLi&Oianl <2014 , on the basis of ultra-precise measure¬ 
ments by the Kepler space telescope, found tiny cyclic variations in 
the period (Light Time Effect - LiTE) of FN Lyr and V894 Cyg that 
are very likely cau sed by a brown dwarf companion. More recently, 
iHaidu et^ j2015ti identified 20 RR Lyrae candidates in binary sys- 
tems in the Galac tic bulge using the same method. Additionally, 
iHaidu et'H] j2015ll proposed tha t about 4 % of all RR Lyraes reside 
in bina ry systems. The results of lLi&Oianl l l2014h and iHaidu et al.l 
( l2015h clearly demonstrated that studying cyclic period variations 
in these overluminous objects is much more efficient in revealing 
binaries than classical methods, i.e. detection of variations in posi¬ 
tion, variations in radial velocities (hereafter RV), or detection of 
eclipses. These results are very encouraging and hold the promise 
for more candidates to be revealed soon. 

In our study we focussed on field RR Lyraes. We 
scanned the literature and the GEOS R R Lyr database] {Groupe 
Europeen d’Observat ions Stellaires, ISoninsegna et al.l 1 20021 : 
iLe Borgne et alj|2007h to choose candidates with cyclic changes 
in their 0~C diagrams which can be interpreted as a manifesta¬ 
tion of an unseen companion. Since there is no review of field stars 
suspected of binarity with reference to various methods currently 
available, we widely discuss particular targets in Sect. As we 
searched for manifestations of the LiTE in 0—C, effects which can 
also affect the appearance of 0—C diagram are discussed (Sect.j^. 
Except for 0—C values from the GEOS database, maxima timings 
from data of various sky surveys and from our own observations 
were determined (Sect. [4). For modelling of the LiTE we used a 
new code introduced in lLiska et al.l (120151) . Results for several tar¬ 
gets are in Sect.^and are generally discussed in Sect.|^ 


2 OVERVIEW OE RR LYRAE STARS SUSPECTED TO 
BE IN BINARY SYSTEMS 

Various types of binary systems with RR Lyrae stars, including 
known candidates, are discussed below. As a reference we pre¬ 
pared a comprehensive list RRLyrBinCan of these objects, includ¬ 
ing false-positives, which is available at the CDS and in a regularly 
updated on-line versior0. 


* http://n'-lyr.irap.omp.eu/dbiT/dbn'-V 1.0_0.php 
^ http://rrlyrbincan.physics.muni.cz/ 


2.1 Eclipsing systems 


This type of binary stars is extremely important, because analysis 
of their light and RV curves allows for the precise determination of 
absolute values of their basic physical parameters. A few RR Lyrae 
stars were proposed to be bound in such systems. Nevertheless, the 
eclipses are in most cases based on unique observations without 
any later confirmation. This could either be some instrumental bias 
or it could be a consequence of a long-period binary with very rare 
eclipses. 

iFitch, Wifoiewski & Johns^ d 19661) reported anomalous be¬ 
haviour of VX Her in one night’s observation data when the star 
was about 0.7 mag fainter than usual in a phase of minimum-light. 
Unfortunately, no repetition of this event has since been reported. 

Similar behav iour was noticed in RW Ari, an RRc type star, by 
IWisniewsO ) ll97 ih . He observed this star during 19 nights in 1966. 
On two nights RW Ari was about 0.6 mag and on another night 
about 0.1 mag fainter than usual in the V-band. He proposed these 
dips in brightness to be primary and secondary minima of an eclips¬ 
ing binary. After subtraction of the average light curve he obtained a 
typical light curve of an eclipsing binary - a very convincing proof 
of his hypothesis (his Fig. 2). According to changes in B —V and 
U—B indices he suggested that the companion is a blu e giant or B- 
type s tar, which is more luminous than RW Ari itself. I^sniewskil 
( Il97ih determined the period of eclipses to be a bit longer than 3 d. 
Another decrease in brightness (0.2 mag, supposed secondary min- 
imum), wh i ch co nfirmed the period of ecl ipses, was re ported by 
IWoodwar^ (Il972h . who analy sed data from Det3 ( 193^. Spectro¬ 
scopic observations made bv lAbt & Wisniewskil ) ll972l) showed a 
very strong difference in RV measured at the same pulsation phase, 
which was interpreted as a consequence of binarity. 

However, lEdwar di ) ll978t) measured RW Ari on 35 nights 
during the seasons of 1976 and 1977 and disproved all formerly 
noted statements about this star. In his observations no indication 
of eclipses was detectable, and the colour of RW Ari was found to 
be very similar to that of other RRc-type stars. He only suggested 
that RW Ari could possibly show Blazko modulation with periods 
of either 38 or 114d. In addition, all following attempts to detect 
eclipses failed. 

It is worth mentioning that reports on eclipses of VX Her 
and RW Ari were based on observations performed by one group 
of observers, at the same time, and using the same equipment 
and methods whereby no comparison star was used (only standar d 
stars, iFitch. Wisniewski & Johnsonll966l : lBookmever et al.ll977h . 
Therefore, the possibility o f a false detection could not be ruled out 
in these two cases (more in lLiska et al.ll201^ . 

V80 UMi, located in the Ursa Minor Dwarf galaxy, is the other 
candidate for an RR Lyrae pu lsating star in an eclipsing system. 
The binary nature, proposed bv lKholopw 1 197lh. is ba sed on pho¬ 
tographic measurements presented in van Agt ( Il968h . Kholopov 
found a possible eclipsing period close to 2 d and a pulsation period 
of about 0.5 d, but both values are dubious, because they are close to 
the integer ratio of one day. In addition the brightness depressions 
were detected only on 10 plates of three consecutive nights. Before 
and after that nothing special was observed. V80 UMi is very faint 
with variations in 6-band of 19.20-20.25 mag, and this was prob¬ 
ably the reason why only one additional dataset was obtained for 
this star dNemec, Wehlau & Mendes de 01iveiralll988h . Their pho¬ 
tographic measurements do not support eclipses. Nevertheless, they 
improved the pulsation period to 0.498746 d and proposed the type 
of the variable to be an anomalous cepheid. 

The systematic search for an RR Lyrae eclipsing system in 
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the Galactic bulge (3256 fundam ental-mode RR L yraes from MA¬ 
CHO project) was unsuccessful ( lRichmondll201 ih . Therefore, the 
only currently known RR Lyrae-like star in an eclipsing binary, 
which was confirmed spectroscopic ally, is OGLE-BLG-RR LYR- 
02792 located in the Galactic bulge dPietrzvhski et al.ll2012h . Un¬ 
fortunately, this low-mass star (0.26 DJI©) with a pulsation period 
of 0.6275 d very likely only behaves similarly to an RR Lyrae pul- 
sator because of its evolution in a close binary systen 0 - the cur¬ 
rent orbital period is 15.24d and the semi-major axis is 32R©. 
The reasons for excludi ng this star from the c lassi cal RR Lyrae 
group were discussed bv IPietrzvhski et aD ( 12012h and lSmolec et al.l 
( l2013h . They found that the bump in the ascending branch of the RV 
curve, Fourier parameter ^31 calculated with theoretical models, 
and a very high rate of period-decrease rules OGL E-BLG-RRLYR- 
02792 out as a regular RR Lyrae star. Therefore IPietrzvhski et al.l 
(l 2012 h defined if as a prototype of a new class of variables, called 
Binary Evolution Pu lsators. Another objec t, OGLE-LMC-RRLYR- 


03541 in the LMC ( 

Soszvnski et alj|2009l). nrobablv also belongs 

to this class dHajdul 

2015. nriv. comm.). However, neither has the 


possibility that the star is a blend of two objects (eclipsing binary 
and RR Lyrae star) been rejected yet. 


2.2 Systems with LiTE 


Binarity of a system with a pulsating star can be unveiled even with¬ 
out eclipses through periodic changes in the times of maximum 
light occurring in the pulsating component, which is caused by the 
orbit of an unseen companion. This is usually examined by using 
0—C diagrams (more in Sect. l3.4t . 

A few RR Lyrae stars were proposed to be in binary systems 
due to LiTE. Probably the most well-known example is TU UMa. 
This star shows sinusoidal propagation in its times of maximum 
light with a period of about 23 years, which is superimposed on 
a parabolic trend caused by ste llar evolution. The idea of bin arity 
of TU UMa was fir st noted by SzeidL Olah & Mizs^ ( 
further anal ysed b y I Saha & White! ( Il99d 3, Kiss et alj ( 

IWade et al.lll999l) . 


198^ and 


and 


Original measurements, current observations available in large 
sky surveys, and observations from photographic plates allowed 
us to perform ver y precise reanalysis of the 0—C of TU UMa 
dUiska et al.ll2015h . It was found that TU UMa very likely com¬ 
prises a pair with a low-mass (minimum value 0.33 9JI©), low- 
luminous star with orbital period of 23.4 years. The binary assump¬ 
tion was supported by the analysis of available RV measurements. 

XZ Dra shows in its 0—C diagram a similar period of cyclic 
variations of about 7200 d (19.7 yr) together with a slow ly increas¬ 
ing pulsation period as analysed bv ijurcsik et al.l ( l2002h . They es¬ 
timated orbital parameters of the possible binary system which al¬ 
lowed them to describe the 0—C diagram and variation in RV mea¬ 
surements. They argued that another mechanism, e.g. a magnetic 
cycle, can play a role. 

Another very int eresting object showin g LiTE-like behaviour 
is BE DoS studied bv IPerekas et al.l l l2004l) . If the binary nature of 
this RRc star is correct, the RR Lyrae component would orbit with 
a period of 8 years around the center of mass together with a black 


^ The secondary, more massive component, is at the beginning of its red- 
gi ant evolutionar y phase . 

^ ISaha & White! dl990t) and IWade et alj dl999ll analysed also systematic 
shifts in RVs. 

^ BE Dor = MACHO* J050918.712-695015.31. 


hole with a mass of about 60 9Jl©. Although they argue against this 
interpretation, and prefer some sort of magnetic cycle instead of 
the binary explanation, the situation with this peculiar object is still 
unclear. _ 

In the last vear. lLi & OianI (1201 4h analysed tiny changes in the 
times of maximum of EN Lyr and V894 Cyg based on data from 
the Kepler telescope. They found that these stars are both accompa¬ 
nied by a low-mass object, probably a brown dwarf or giant planet, 
with orbital periods of 795 d (FN Lyr) and 1084 d (V894 Cyg). It is 
very likely that mor e such binaries will be unveiled i n precise space 
measurements (e.g. lGuggenberger & Steixneill2014l) . 

Very recently, twenty binary candidates w ere identified among 
OGLE bulge fundamental-mode RR Lyraes bv iHaidu et al.l ( l2015h . 
They started analysis from 1952 of well-observed stars with base¬ 
lines of up to 17 years. From these results, they estimate that > 4 % 
of all RR Lyrae stars belong to binary systems. 

Analysis of the long-term per i od behaviour of AV Peg, RS 
Boo, RR Leo dOlah & Szeidlljl978l. iSzeidl. Olah & MizseJll^ 
and RZ Get dLe Borgne et alj 2007h also provid ed some indica¬ 
tion o f the binarity of these objects. In addition. iLe Borgne et al.l 
( l2007h gave several other somewhat tentative examples with some 
peculiarity in their 0—C diagrams, which can possibly be a conse¬ 
quence of binarity (e.g. SX Aqr, X Ari, etc). Long-term monitor¬ 
ing, lasting a few decades, facilitated the identification of several 
RR Lyrae candid ates in binary systems in globular clusters, e.g. 
M3, M 5 or M53 ( ISzatmarvlll9mlSzeidl et alj|201ll:ljurcsik et al.l 

ImI. 


2.3 Spectroscopic binaries 

Spectroscopic binaries, discovered through variations in RV, com¬ 
prise a large fraction of known binary systems among ordinary 
stars. Nevertheless, this statement, at least partially, does not apply 
to RR Lyrae stars, in which the situation is complicated due to pul¬ 
sations. In addition, in detached systems, the pulsation component 
probably would dominate in the spectrum (reasons were mentioned 
in Sect.O. Thus, only long-term, high-precision spectroscopic ob¬ 
servations, and application of spectral disentangling, could bring 
long-desired success. 

So far the only known spectroscopic binary with an RR Lyrae 
comp onent is OGLE-BLG-RRLYR-02792 dPietrzvnski et al.l 
l2012h . which is also an eclipsing binary (see Sect. Ol. This star is 
a double-line spectroscopic binary with the main component only 
mimicking an RR Lyrae type star. Concerning galactic field stars, 
TU UMa is one of the most promising targets, in whi ch possible 
RV c hanges should be observable during one decade dUiska et al.l 
20lj). Variation i n RV was also analysed in XZ Dra (Sect. 12.21 
Jurcsik et alj2003) . 

Several candidates for spectroscopic binaries were identified 
based on differences in systemic RV from different studies which 
achi eve up to several tens of k m s^*. The list of suspected systems 
from iFernlev & BamesI dl997h contains Cl And, DM Cyg, BK Dra, 
XX Hya, ST Leo, CN Lyr, and TU UMa, and is co mplern ented 
by TY Aps, BX Dra, and BX Leo from I Solano et al.l dl997h . DM 
Cyg and TY Aps show Blazko modulation, which could produce a 
shift in RVs, and BX D ra was re-classified as an eclipsing binary 
dAgerer & Dahmiri995h . 

2.4 Visual binaries 

Although the great distances to RR Lyrae stars make visual pairs 
hardly detectable, this type of binaries should be mentioned at least 
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for completeness. The upper limit for the separation of gravitation¬ 
ally hound objects is prohahl y not known. Never theless, from a 
catalogue of 3139 orbits from iMalkov et ^ ll 2012 h . we found that 
99.3 % stars in 2869 binaries have semi-major axes shorter than 
1000 au. For example the prototype RR Lyrae itself (the closest RR 
Lyrae star, parallax tz = 3.46 mas. Ivan Leeuwenll2007h as a hypo¬ 
thetical system with semi-major axis of 1000 au would have a com¬ 
ponent with a maximum angular distance of 3.46". This illustrates 
the difficulty in finding visual systems among RR Lyraes, because 
possible visual systems would be much closer and changes in posi¬ 
tion angle undetectable on reasonable time scales. The chance for 
discovery of such a system without adaptive optics or interferome¬ 
try is negligible. 


cepted ephemeris. Subsequently the time difference 0—C is plotted 
against time or a cycle number N. 

In an 0—C plot it is very easy to investigate long-term sys¬ 
tematic trends, which can be e.g. monotonically linear (constant 
period) or parabolic (period change with a constant rate), periodic, 
as well as sudden, more or less ir regular (a review of 0 —C dia¬ 
gram shapes can be found e.g. in iLiska & Skarka1l2015h . In the 
case of RR Lyrae stars, the behaviour of their periods can be very 
complex with many possible reasons combined, sometimes with an 
unclear interpretation. Therefore, an RR Lyrae star showing LiTE, 
without another independent confirmation, should only be consid¬ 
ered to be a “candidate” binary system. Very illustrative examples 
of many 0—C diagrams c an be found e.g. inijurcsik et al.l ( 1200 iL 
l 2012 h : ISzeidl et aLfeoilhJLe Borgne erS](l2007h . 


2.5 Other types of binarity 


In the second half of the twentieth century, AR Her and 
BB Vir were found to be sign ificantly hotte r at minimum 
light than other RRab vari ab les iPrestonI 1 9591: Stur^ 196^ 
iFitch. Wisniewski & Johns^ 1 196a Bookmever et a ll Il977h . 
Later, this discre p ancy was discusse d and explained by 
iKinman & Carrettal ( Il992h and iFernlevI ( Il993h as a possible 
consequence of binarity. They considered the second component 
to be a blue, hot star. 

Both BB Vir and AR Her are known as Blazko stars, but the 
modu l ation nature of BB Vir is somewhat questionable dSkar 
l2014h . iKinman & Carretti (Il992h analysed photometric observa¬ 
tions of BB Vir in detail using (B-V) colour and also light curve 
characteristics. They found that BB Vir has a lower (B-V) than 
other RRab stars, and that the amplitude of its light curve is smaller 
than it should be for the given colour. After their correction of the 
BB Vir light (assumed to be similar as for other RRab stars), they 
suggested the unresolved component to be a hot HB star with sim¬ 
ilar physical characteristics as BB Vir itself, which lies on the blue 
edge o f the instability strip. Based on ultraviolet spectra, iFemlevI 
(Il993h more closely specified the temperature of the supposed com¬ 
panion to 7900 K. He also suggested an alternative explanation of 
BB Vir to be an RRc st ar with an unusually long period. 

iBorkows^ (Il980l) explained the modulation of AR Her to be 
the consequence of non-linear superposition of the fundamental 
(P] = 0.4700 d) and the second or third overtone (P 2 = 0.233 d). 
IKinman & Carretti d 19921) speculated that, if the binary nature of 
AR Her is correct, the contribution of the second overtone to its 
modulation could actually be explained by the presence of an un¬ 
resolved RRc star with a 0.233-day period r ather than the con¬ 
sequen ce o f beating betw een radial overtones. IKinman & Carrett^ 
dl992l) and iFemie^ d 1 993l) nicely showed how expected photomet¬ 
ric and spectroscopic characteristics of RR Lyrae stars would be 
affected by a possible companion. Despite all the discussed theo¬ 
ries, a convincing proof of the binarity of BB Vir and AR Her is 
still missing. 


3 PERIOD CHANGES OF RR LYRAE STARS IN 
O C DIAGRAM 

Since the beginning of modern variable star research, the very sim¬ 
ple and illustrative 0—C diagram has been used for the description 
of the period evolution of stars. Time of a given phase (usually max¬ 
imum or minimum of the light changes), which is observed (O), is 
compared with predicted time calculated (C) on the basis of an ac- 


3.1 Evolutionary effects 

Already at the beginning of the last century, lEddingtc^ dl918l) 
pointed out that pulsating stars should undergo evolutionary pe¬ 
riod changes due to changes in their density. These changes are 
very slow, therefore observations spread out over many decades are 
needed. Eortunately, RR Lyraes have been observed for more than 
a century in some particular cases. 

In RR Lyrae research, it is usual to express the period-change 
rates in j3 = Ppyjs = dPpyjs/df in [dMyr^*] or [msd^*], and 
a = j6/Fpuis = (l/fpuls)^puls [Myr^*] parameters, for definition 
of Ppuis = 2 a 3 /Ppuis, where 03 is a third coefficient in the parabolic 
ephemeris of maxima timings (e.g.. iLe Borgne et al. l l2007lFi De¬ 
pending on their masses, chemical composition and evolutionary 
stage, RR Lyrae variables can exhibit both a period decrease and a 
period increase during their evolution on the HB (see e. g. the theo¬ 
retical works of lSweigarlll987l : lLee & Demarauelil990l) . When the 
evolution is blueward, the period of a star slowly decreases, while 
the evolution from blue to red implies a period-lengthening. Since 
the redward evolution is more rapid than towards the blue, the rate 
of change of the period is higher for period-lengthening than for 
period shortening. _ _ 

According to iLee & Demarguel ( Il99(ll) . the majority of RR 
Lyrae stars are in the early stage of their HB evol ution. Analysing a 
large s ample of stars such as in globular clusters. fLee & Demargud 
d 19901) found that the average value of j 6 should be around zero 
or ver y slightly positive. This assumption was confirmed by iLed 
( Il99ll) . who compared synthetic with the observed HB and found 
that jS depends on the type of HB . An other confirmation came with 
studies ljurcsik et ^ ( I200ll . l2012l) and lSzeidl et alj ( l201l[) . that also 
investigated period changes of RR Lyraes in globular clusters. 

Nevertheles s, _ in all these studies, as well as e.g. in 

iLe Borgne et alj | |2007t) . who analysed the 0—C of field stars, some 
peculiar cases with a very high positive or negative j 6 were re¬ 
ported. Stars with a large j3 > 0 are considered to undergo the 
final episodes of their HB life, while RR Lyraes with the most 
rapid period-decrease are s upposed to be in their pr e-HB stadium 
dSweigart & Renzitiill 19791 : Isilva Aguirre et alj|2008l) . However, in 
several cases with the largest /3, such high period changes can¬ 
not be explained by evolution theory. We suppose that for some of 
these stars the discrepancy between observation and theory could 
be caused by misinterpretation of their 0—C diagrams. In the case 


® Several authors fe.g. iHaidu et al.|[2015l) used different definition of /l- 
parameter - their /Ifhcir = <^3 and then jSthcii- = 1 /^PpuisP- 
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of a potential double star with a long orbital period, the correspond¬ 
ing 0—C could mimic a parabola when only a short time span of 
observation is available. 

In a few objects the shape of the 0—C diagram is more com¬ 
plicated than only a simple parabola. A cubic ephemeris with sig¬ 
nificant third-order terms was also someti mes calculated for thes e 
objects (e.g., RS Boo, SU Dra, RW Cnc, iLe Borgn e et al.ll2 007h. 
This a pproach is meaningless as was shown by Mikulasek & Zeida 
( l2013h . Evolutionary changes described using the power law, in 
which dPpuis/d/ ~ , where q is the deceleration parameter, the 

third-order term is negligible even in objects with extreme period 
shortening/lengthening. This implies more complex behaviour than 
simple evolutionary effects are present in these objects (e.g. bina- 
rity). 


3.2 Blazko effect 

Almost a half of all RRab stars l ljurcsik at alJl200^ : iBenko et al.l 
l2010h shows more or less periodic changes of the shape of thei r 
light curves, which is known as the Blazko effect Hi ^ Il907h . 
These changes manifest as changes in amplitude and in time 
of maximum light. Therefore, these changes can be reproduced 
and analysed using an 0—C diagram. An 0—C diagram of a 
Blazko star shows a periodic pattern with an amplitude of maxi¬ 
mally about 0.2 d (the c urrently known record hol der is V445 Lyr 
with A(0—C)~ 0.15 d. iGuggenberger et'^l2012h . Periods of the 
Blazko effect differ by several orders - from 5 days to more than 
5 years dSzczvgiel & Fabrvckvll200^ . In fact, the upper limit for 
the length of the Blazko cycle is not known, and the greatest values 
detected for the period are governed by the length of the time-base. 
Thus the Blazko effect can be easily misclassified with LiTE. 


3.3 Erratic and other changes 

Except for periodic changes caused by the Blazko effect, RR 
Lyraes undergo very sudde n, random changes, which affect 
mainl y Blazko stars (see e.g. ljurcsik et alj|2012 : iLe Borgne et al.l 
l2007h . T hese abrupt, significa nt changes o bserved, for example, in 
XZ Cyg l lBezdenezhnvil 1 988h or RR Gem jSodor. Szeidl & JurcsikI 
l2007h . are usually associated with instabilities during the final 
phase of the helium-burning phase or a s a consequence of the mix - 
ing events at the convective-core edge dSweigart & Renzinilll979h . 
Other possible explanations deal with changes in the gradient of the 
helium compositi on in the re gions below the hydrogen and helium 
convection zones l lCoxlll998h . 

Additional, not necessarily erratic changes, can be connected 
with hydromagnetic effects. They are probably the only way to ex¬ 
plain the problem with high-amplitude cyclic variations detected 
in non-Blazko stars which appears like peculiar LiTE caused by a 
very massive com ponent (a possible bl ack hole with several tens of 
2)?0, e.g. BE Dor. lDerekas et al.ll2004h . 


3.4 Light Time Effect 

Probably the last possible way to affect the 0—C values of RR 
Lyrae stars, is through the effect of orbital motion in a binary sys¬ 
tem with appropriate parameters (large semi-major axis and incli¬ 
nation). The orbital motion of an RR Lyrae component will pro¬ 
duce detectable changes of its times of maximum light. This is a 
consequence of the finite speed of light, and it is called the Light 
Time Effect (LiTE), as already mentioned. When a pulsating star 


Table 1. Number of new times of maxima determined from individual 
projects (HIP - Hipparcos), from our observations, and GEOS database. 


Star 

ASAS-3 

CRTS 

DASCH 

HIP NSVS OMC 

SWASP Own GEOS 

Interval 

RS Boo 

0 

0 

49 

3 

4 

0 

32 

0 

448 

1900.0-2013.5 

RU CVn 

0 

5 

24 

0 

6 

0 

43 

0 

114 

I902.0-20I3.3 

RZCet 

15 

8 

0 

3 

9 

4 

0 

0 

100 

1930.0-2013.8 

S Com 

10 

10 

34 

3 

11 

0 

0 

0 

134 

1901.2-2013.4 

SU Dra 

0 

0 

0 

5 

5 

0 

0 

0 

235 

1904.2-2013.5 

VX Her 

13 

4 

0 

4 

9 

0 

0 

6 

243 

1916.2-2014.6 

SS Leo 

13 

2 

0 

2 

3 

0 

0 

3 

105 

1901.3-2014.4 

AVPeg 

9 

0 

0 

3 

4 

0 

0 

0 

502 

1931.7-2013.6 

ATSer 

14 

7 

0 

3 

6 

1 

0 

3 

50 

1900.9-2014.5 

RV UMa 

0 

0 

53 

5 

14 

12 

0 

0 

252 

1897.3-2013.3 

BB Vir 

5 

7 

38 

2 

5 

6 

0 

0 

46 

1902.1-2013.4 


is in the closest part of its orbit to the observer, the time needed 
to reach the observer is logically shorter than when it is in other 
parts of the orbit. As the star revolves around the center of mass, 
the corresponding 0—C plot will change periodically. The shape 
of these changes depends on the orientation and eccentricity of the 
orbit. The amplitude of such an 0—C diagram directly provides in¬ 
formation about the projected size of the orbit, and the period of 
these changes correlates with the orbital period. These parameters 
allow us to estimate the total mass of the system. A few systems 
suspected of LiTE were briefly discussed in Sect. 12.21 


4 DATA SOURCES AND TARGET SELECTION 


Since the GEOS RR Lyrae database dLe Borgne et al.ll2007h is the 
most extensive archive containing times of maxima of RR Lyrae 
stars, it was used as the main source of data for this work. Targets 
from literature suspected of LiTE were supplemented by objects se¬ 
lected only on the basis of visual inspection of their 0—C diagrams 
from the GEOS web page. This means that only stars with suspect 
cyclic variations were selected (stars showing one or more cycles 
or almost one cycle dissimilar to a parabola). In addition, only stars 
with more than 10 maxima timings passed for further analysis from 
the entire GEOS database. After a closer inspection of the remain¬ 
ing 0—C diagrams only eleven RRab stars, which are discussed in 
Sect.|5] were fully analysed. Erom the preceding discussion it is ev¬ 
ident that we focused on high amplitude and long-term variations 
in 0—C diagrams with good time-coverage. Therefore, our study 
is strongly affected by selection effects. 

We did not use data marked as ‘pr. com.’ from the GEOS 
database (the numbers of used maxima from the GEOS database 
are in the Table [TJ. To maximally extend the GEO S sample w e 
re-anal yseciZi data from sky surveys - Hipparcos (|ESA| 1997), 


NSVS dWt^niak etal .112004 . ASAS-3 jPoimanskil 119971 

i 


2002 ), 


and g ained new ma xima from CRT^ , DASCH jOrindl av et al.l 
l2009l) . SunerWA SP jPollacco et al.ll200d:lButters et al.ll2010h. and 
OMC ^ (see e.g lAlfonso -Garzon. Domingo & Mas-Hesse 20 id : 

3 


I Zeida^ Doming J 201 ih . New maxima for the studied systems 
were determined from sparse data using the template fitting method 


^ From the reasons discussed in lLiska et alJi2Ql3) . 
^ http://crts.caltech.edu/ 

^ http://sdc.cab.inta-csic.es/omc/index.jsp 
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described in iLiska et ^ l l2015h . The dataset with the best qual¬ 
ity (good accuracy, phase coverage) was selected for each star for 
the template fit. Maxima from SuperWASP measurements (good 
time-resolution) were calculated using polynomial fitting. The de¬ 
termined times, including uncertainties, obtained directly from the 
least-squares method are available through CDS. The number of 
used maxima from particular surveys can be found in Table [T] 

New photometric measurements in BVRI bands were per¬ 
formed for VX Her at Masaryk University Observatory (MUO), 
Brno, Czech Republic during 13 nights (April-August 2014) with 
a 24-inch (0.6-m) Newtonian telescope and a G2-0402 CCD cam¬ 
era. Other measurements in phases near maximum brightness for 
the stars SS Leo and AT Ser were also obtained in Brno with a 
small 1- inch telescope in the g regn-band (description of the tele¬ 
scope in lLiska & Liskovil2014 1. 

All the maxima used from the GEOS database, with corre¬ 
sponding references, maxima determined from sky surveys mea¬ 
surements, and from our observations, are available via CDS. For 
each star the interval covered by these measurements is given in 
Tabled 


5 RESULTS ON STARS THAT ARE SUSPECT 

The process of LiTE fitting described in lUiska et H] ( 1201 Sh was 
applied on our selected target stars whose identification and basic 
light ephemeris are in Table[2 Each of these stars shows significant 
period changes. Several of them were previously reported as possi¬ 
ble members of binary stars, but with undetermined orbital periods, 
or only very roughly estimated periods. 

Some of the 0—C changes can be interpreted directly as large 
amplitude cyclic changes. Such dependences were fitted as a whole 
and are described by our model 1. In several cases the orbital pe¬ 
riod is similar to the length of the time-base. Therefore model 1 
should be considered only as a possible solution, in which the ob¬ 
served changes may not be cyclic at all. Some of the stars’ 0—C di¬ 
agrams can be interpreted as a low amplitude cyclic change super¬ 
imposed on a secular period change. This situation is described by 
our model 2. However, in some cases both interpretations could be 
possible and, therefore, both models were computed. Our results 
with the parameters describing the proposed orbits are in Table [3] 
The columns are described in the notes below the table. Except for 
orbital parameters this table also contains the period-change rate 
(if present), goodness of fit indicator (Xr) ™d the number of used 
maxima timings. For calculation of minimum masses of compo¬ 
nents we assumed an RR Lyrae mass of 0.6DR©. 

In agreement with expectations, it seems that companions 
of the majority of our sample stars are low-mass stars. However, 
in RZ Get the companion is probably a white dwarf or neutron star, 
and in AT Ser it is probably a black hole. The studied stars are 
discussed in detail in the next sections. 

Available literature contains a relatively rich sample of origi¬ 
nal RV measurements (Table 0 which could possibly be used for 
verification of our results, because our models allow the computa - 
tion of systemic variations in RV (for details see lUiska et allTOfl) . 
Unfortunately they are in most cases of an insufficient number, time 
span, or quality for such usage. Systemic veloci ties for our tar¬ 
gets fr om iLavdenI (ll994h . lFernlev & BamesI h997h or lSolano et al.1 
( Il997l) . which are typically based on 2 or 3 unpublished measure¬ 
ments, are of insufficient quality. During our study we found that 
among our targets some of th e RV measureme nts may not be reli¬ 
able - mainly old values from lAbtl d 19701 . [197^ can differ by more 


Table 4. Sources of radial velocity measurements for individual RR Lyrae 
stars. 


Stai- 

Ref 

Star 

Ref 

RS Boo 

2, 8, 14 

SS Leo 

1,3, 5, 7 

RUCVn 

2 

AVPeg 

2, 9, 14 

RZ Get 

4 

AT Ser 

2, 6,7 

S Com 

2,7 

RVUMa 

2, 6, 12 

SUDra 

1,9, 10, 11, 13 

BB Vir 

2 

VX Her 

2, 14 




J 1 9.5(11. 5 -iFemlev et al 

h99(]tl. 6 -iFernlev. Skillen & Burkil 

1993h. 

7- Hawlev & BamesI Jl985l). 8 -IJones. Camev & LathamI Jl988 

1. 9- 

iLiu & Janed h989ll. lO-lOke. Giver & Searle] h962[). 11- 

Preston! 

196.5h. 

12-Preston & Sr)inradHl967h. 13- Varsavskvl ll96Cll. 14- 

Woollev & Alvl 


than 30 km s * compared to other dataset This shows the need 
for new accurate spectroscopy. 


5.1 RS Boo 


RS Boo is one of the best studied RR Lyraes. Except for long¬ 
term period changes, it shows modulatio n (Blazko effect) with a 
period of about 533 d l lOosterho fllll946l) and a possible shorter 
period of between 58-62d (lKanvdll980h . More recent analy- 
ses showed multiple B lazko modulation with p eriods of 532 .481 d 
dLe Borgne et al.|[2012h and 41.3 d and 62.5 d ( Skar^|2014h. The 
idea that RS Boo could be a binary comes from Kanvdl ( 198(d) who 


investigated its 0—C diagram. After subtraction of a parabola he 
proposed an orbital period of about 70 years based on residuals. 
I Kanvdl dl986h also found that the brightening during the 62-day cy¬ 
cle is lower in blue than in yellow, contrary to the long 530-day 
cycle. He guess that it could be due to a cooler companion. 

The shape of the 0—C diagram of RS Boo (significantly ex¬ 
tended with values based on the DASCH project) is apparently 
asymmetric (Fig.Q} allowing direct interpretation as a consequence 
of LiTE (our model 1) and as a secular period lengthening with 
superimposed LiTE (model 2, Fig.[T]l. Model 1 gives an orbital pe¬ 
riod of 147 yr and a minimum mass of the companion of 0.264 DJI©, 
while model 2 gives an orbital period of 77 yr and a minimum mass 
of the second component of 0.097 DJI©. Although the shape of the 
0—C also suggests non-symmetrical period lengthening^ the in¬ 
dication that the dependence should go makes this explanation less 
plausible. Because both models have the same quality of the fit 
(X^ = 1.082) and both secondary component masses are possible, it 
is very difficult to choose the correct solution. Only further precise 
photometric and spectroscopic measurements could help. 

After subtraction of both models we searched for additional 
periodicity in residuals. Both models give almost the same resid¬ 
uals, therefore the same periods were found in both cases. The 
most significant peak in the frequency spectrum (top panel of the 


iLaydgij (1994 ) anal ysed a large sample of mean RVs from 
IPavne-Gaposchkir il954h for RR Lyraes and estimated the uncertainty of 
these measurements as 35kms“*. In many cas es original values for these 
m ean RV measurements w ere published later in lAbd dl970l[l973h . 

** iLe Borgne et al ]i2^ performed a third-order least-squares fit. 
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Table 2. Identification and basic light characteristics of studied stars (taken from The International Variable Star Index - VSX. IWatson et alj2QQd) . 


Star 

RA 

DEC 

Other designation 

Mo [HID] 

Tpuis [d] 

Mag range (V) 

RS Boo 

14:33:33.21 

+31:45:16.6 

AN 14.1907 = HIP 71186 

2441770.49 

0.37733896 

9.69-10.84 

RUCVn 

13:59:33.28 

+31:39:04.3 

AN 6.1914 

2434483.467 

0.5732449 

11.36-12.48 

RZ Cet 

02:28:32.44 

-08:21:30.1 

AN 36.1929 = HIP 11517 

2433906.8920 

0.5105981 

11.245-12.240 

S Com 

12:32:45.63 

+27:01:45.4 

AN 83.1910 = HIP 61225 

2440654.6410 

0.58658722 

10.88-12.12 

SU Dra 

11:37:56.61 

+67:19:47.0 

AN 43.1907 = HIP 56734 

2443902.0467 

0.66042001 

9.18-10.27 

VXHer 

16:30:40.73 

+ 18:22:00.2 

AN 30.1917 = HIP 80853 

2453470.2068 

0.4553595 

9.913-11.177 

SS Leo 

11:33:54.49 

-00:02:00.0 

AN 69.1919 = HIP 56409 

2441781.4090 

0.6263351 

10.420-11.584 

AVPeg 

21:52:02.79 

+22:34:29.4 

AN 90.1931 =HIP 107935 

2443790.3160 

0.3903814 

9.93-10.99 

ATSer 

15:55:40.44 

+07:59:18.7 

AN 46.1935 = HIP 77997 

2453467.0077 

0.7465617 

11.013-11.890 

RVUMa 

13:33:18.09 

+53:59:14.6 

AN 139.1907 = HIP 66122 

2445075.511 

0.46806 

9.81-11.3 

BB Vir 

13:51:40.78 

+06:25:51.4 

AN 23.1935 = HIP 67653 

2439613.7760 

0.4711069 

10.732-11.450 


Table 3. Pulsation and orbital parameters for studied systems with probable LITE*. 


Star 

Model 

^ puls 

Mq 

10-'‘f’p„i, 


^ 3rbit 

To 

e 

ft) 

A 

sin/ 

/(ffl!) 

®t2^niin 

Kl 

Zr 

^max 



[d] 

[HJD] 

[dd-‘] 

[dMyr-l] 

[yr] 

[HJD] 


[°] 

[light day] 

[au] 

[ffltol 

[fflto] 

[kms“* ] 



RS Boo® 

1 

0.37733735t]) 

2448500.31931)) 

- 

- 

1471)5 

24255001)555 

0.1501“ 

2201)’ 

0.04681)) 

Q 1 -i-1.3 
^•‘-0.8 

0.02461)* 

0.26411) 

1 . 6601 !! 

1.082(61) 

536 

RS Boo® 

2 

0.3773376891)1 

2435747 . 03461)5 

+ 12.811)5 

+0.04681)* 

16.1^21 

24553501)* 

0.551)) 

2131)1 

0.01291)'* 

2.231)® 

0.001881)5 

0.0971)’ 

1.031’® 

1.082(62) 

536 

RU CVn 

1 

0.5732503141^® 

2455343.46121)) 

- 

- 

101.61);) 

244850015)55 

O. 22 I)* 

3331)5 

0.04951)) 

8.615) 

O.O 61 I* 

O. 392 I)! 

2.581)5 

1.03(10) 

192 

RZCet®’ 

2 

0.510612441)5 

2441276.62481®) 

-180.8l);5 

-0.6601)5 

'^•^-4.0 

243153015“ 

0.351)® 

105.21®;! 

0.082115 

14.111:’ 

0.4911! 

1.151)) 

5.941’® 

1.13(12) 

139 

S Com 

1 

0.5865890921* 

2453796.555381” 

- 

- 

105.815) 

24422001)555 

0.1 2 I)* 

1351)) 

0.03291)* 

5.701)5 

0.0166115 

0.22411)5 

1.6171* 

1.05(10) 

202 

SCom 

2 

0.58658912415’ 

2435949.00671)) 


- 0 . 0261)5 

90.ll);| 

24427001)55" 

0.1931)® 

1341®’ 

0.02281)) 

3.961)) 

0.00761)) 

0.1651)) 

1.3311’ 

1.03(10) 

202 

SU Dra 

2 

0 . 660421241)2 

2443902.05491)) 

+ 16.61);) 

+0.06061)) 

70.71);) 

24321001)555 

0.241)’ 

1091* 

0.01131)) 

1.961)5 

0.001501)5 

0.089115 

0.851)® 

1.118(91) 

245 

VX Her®- 

2 

0.455366782l)( 

2438911.56831)4 

- 42 . 71 I 55 

-0.15601)) 

83.01);5 

242210011555 

0 . 8 III* 

1341!® 

0.007871)) 

1.3611® 

0.000371)* 

0.05401®* 

0.8411) 

1.095(86) 

279 

SS Leo 

1 

0.626341331)5 

2441781.391411) 

- 

- 

110.71)) 

24532401))5 

0.4751)5 

141115 

0.05171)5 

8.951)5 

0.05851!) 

0.3841)5 

2 . 737 I 55 

1.09(13) 

128 

AVPeg 

2 

0.3903752821)) 

2441552.288771)5 

+45.961)) 

+ 0 . 16791)5 

47.7l);| 

24558801)“ 

0.441)* 

86115 

0.005911)5 

1.0241)5 

0.000471)® 

0.058915! 

0.711)* 

1.010(63) 

518 

AT Ser 

1 

0.7465600311) 

2436093.26091)5 

- 

- 

85.611) 

24494401)55 

0.4561)5 

217.71);® 

0.11571* 

20.015-5 

1.101)5 

1.901)5 

7.831)5 

1.12(16) 

84 

RV UMa® 

2 

0.4680628501)) 

2438459.9318115 

+4.7015) 

+0.01721)) 

66.9315) 

24495201)55 

0.4041* 

3171* 

0.026571)5 

4.60115 

0.02171)5 

0.2511)1 

2.241)1 

1.105(78) 

336 

BB Vir®’ 

2 

0.471099581), 

2436232.24261)5 

+36.311;) 

+0.132515'* 

92.71);’ 

244690015)55 

0.241]* 

2961)5 

0.02481J) 

4.30l’l 

0.00921)) 

0.1771)) 

1.421)5 

1.11(14) 

109 


Notes, h) Columns contain following parameters: Star - name of the star in the GCVS, Model - type of calculated model (1 = only LITE, 2 = LITE + parabola), Ppuis - pulsation period, Mq - zero epoch 
of pulsations, Ppuis = ^ - relative rate of changes of pulsation period, Porbii “ orbital period, Tq - time of periastron passage, e - numerical eccentricity, CO - argument of periastron, A - sin / in light days 
(semi-amplitude of LiTE can be calculated as An-pg = A \J\—e^ cos- ft)), sin/ - projection of semi-major axis of primary component according to the inclination of the orbit /, fipyi) - mass 

function, VJ12 min “ lowest mass of the second component, the value was calculated for inclination angle i = 90° and adopted mass of primary fDlj = 0.6 OJt©, Ki - semi-amplitude of RV changes primary 
component (RR Lyrae star), - normalised value of where “ g) for number of used measurements A'max and number of free (fitted) parameters g (only LiTE g = l, LiTE -l- parabola 

g = 8), A'max - number of used maxima timings, - Blazko effect is known in the star, more information in the text, - the star is suspected from Blazko effect, more information in the text. 


Fig.E) relates to a Blazko period of Bbl = 535(4) d (for all avail¬ 
able data), which is within the error margin the same period as 
from previous studies (e.g. lKanvdll984 ke Borgne et alJIlOOl 
l2012h . and PgL = 550.5(1.0) d when only more recent photoelec¬ 
tric and CCD measurements were used (data phased with this pe¬ 
riod are in Fig.[3l. After prewhitening with this peak, no other sig¬ 
nificant peak was revealed taking all measurements into account. 
On the other hand, when frequency spectra of CCD and photo¬ 
electric data were investigated, three additional peaks were iden¬ 
tified. The first of them, with unclear interpretation, corresponds 
to a period of 2400 d. The other two peaks, corresponding to pe¬ 


riods P^2 = 39.5(7) d and = 60.0(1) d (two middle panels in 
Fig. ID, correlate well with the secondary modulation periods from 
literature. Therefore we argue that the secondary modulation com¬ 
ponents are real and need to be refined by additional observations. 

Our first model of LiTE suggests RV variations with a semi¬ 
amplitude of Ki ~ 1.7 km s^* (or only Ikms^* from model 2) 
which is too small for reliable confirmation with the available RV 
measurements (see Table ID. It is becaus e the only measurements 
with a sufficient accuracy, the ones of Ijones. Carney & Lathanil 
(19881) . span only 108 days in various phases of the 550-day cy¬ 
cle, and therefore they are useless for our purposes. At least these 
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RS Boo model 2 



0.12 
0.1 
0.08 
0.06 
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20000 30000 40000 50000 60000 

HJD-2400000 [d] 


RS Boo model 2 - parabola subtracted 



10000 20000 30000 40000 50000 60000 

HJD-2400000 [d] 


Figure 1. O—C diagram of RS Boo together with our model 2 (solid 
line, top panel) and variation after subtracting the pai'abolic trend (bottom 
panel). Black circles and blue stars display maxima adopted from the GEOS 
database and new maxima determined in this work (from sky surveys mea¬ 
surements and from our observations), respectively. 


data allowed us to reject the possibility that the 550-day cycle could 
be a consequence of an additional LITE, because such LITE would 
produce periodic variations in the systemic velocity with a semi¬ 
amplitude of 20. 9 km s~' which, after subtracti on from the val¬ 
ues adopted from ljones. Carney & Lath^ l ll988l) . results in an in¬ 
crease in the scatter of the phased RV curve. 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 




Frequency [c d'-l] 


Figure 2. Frequency spectra of photoelectric and CCD maxima for RS Boo. 
Steps in prewhitening with labelled frequencies go from the top to the bot¬ 
tom panel, where residuals are shown. 
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5.2 RU CVn 

Possible long-term cyclic changes in the pul s ation p eriod of prob¬ 
ably 60000 yQ was discussed by iHusad ( I 2 OO 3 I) . and due to 
the complicated shape of the 0—C diagram was mentioned in 
iLe Borgne et alj ( |2007|) . The dependence in Fig. |4] clearly shows 
a cyclic variation. Thus only model 1 was applied, which gives an 
orbital period of about 102 yr and a minimum mass of the compan¬ 
ion of 0.39 tPl©. However, orbital parameters based on our model 
should be considered only as preliminary, because only one cycle 
has been completed up to now. 


This value should probably be in days not in years. 


Figure 3. Residual 0—C values of RS Boo after subtracting model 1 phased 
with the found Blazko period of 550.5 d. Vertical scatter is caused by addi¬ 
tional modulation components. 

R U CVn has only 2 single RV measurements available in lit¬ 
erature l lAbtll973h . Both values were obtained with a nine-year dif¬ 
ference at nearly the same pulsation phase of 0.33, but they differ 
significantly (-1-21.7, —58.9) kms^*. The big difference in RVs can¬ 
not be explained by the difference in systemic velocity based on our 
LiTE model, because changes due to two different phases of the or¬ 
bit should only be about 1 kms^*. In addition, the semi-amplitude 
of RV changes K\ is only 2.6kms^*. One of the values is proba¬ 
bly incorrect or values are highly scattered (this is a more plausible 
explanation), see Sect.[3 
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RU CVn model 1 



10000 20000 30000 40000 50000 60000 


HJD-2400000 [d] 

Figure 4. O—C diagram of RU CVn together with our model 1. Symbols 
are the same as in Fig.[T] 


5.3 RZ Cet 

Period variations of RZ Cet in its 0—C di agram are well 

S ironou nced (Fig. O and were studied e.g. by iLe Borgne et al.l 
2007h . They fitted the 0—C dependence with a parabola (Ppuis = 
— 157.9(8.8) X 10^" dd^*) and subsequently found oscillations in 
residuals with a period of 12500 d superimposed on the parabolic 
trend. They mentioned that this oscillation could be the conse¬ 
quence of LiTE. 

We tested the proposed scenario using a similar approach - 
simultaneous fitting of a parabola and LiTE and found Ppuis = 
-180.8+^'^ X 10^" dd^*, and a LiTE period of 75yr (27600d) 
which is more than two times more than the value from 
iLe Borgiie et alJ (l2007h . Nevertheless, the time-baseline has a sim¬ 
ilar value (84 yr) and thus the orbital period and other parameters 
(such as companion mass 97l2.min ”^1.15 DJte) are very prelimi¬ 
nary. 

However, variations in the 0—C diagram allow for another, 
less probable scenario: that the 0—C dependence itself could be 
interpreted as a LiTE (without parabola, model 1) with incom¬ 
plete cycle. The insufficient time span caused the modelling to be 
unstable and the LiTE model parameters had to be manually ad¬ 
justed. This ’solution’ gives an orbital period of at least 86 years 
with an amplitude suggesting a high-mass companion with a min¬ 
imum mass of 7.5 DJlp,. According to the discussion in Sect.[T]this 
should be a black hole. Since the proposed orbital period is more 
than 80 years, the system would be well detached and therefore no 
mass transfer would be possible, and the black hole would manifest 
itself only via gravitational effects on the RR Lyrae component. Al¬ 
though the presence of a black hole could be surprising, RZ Cet 
would not be the first object among RR Lyrae type stars with a 
poss ible black hole (see the discussion about BE Dor in Sect. 12.21 . 

iKovacd ( l2005h identified RZ Cet as a modulated s tar, but 
the pe riod of the proposed Blazko effect is not known. Iskar 
1 I 2 OI 4 I 1 analysed data of RZ Cet from the ASAS survey and found 
no indications of modulation. It is possible that high-amplitude 
0—C changes could be attributed to the Blazko effect with an ex¬ 
tremely long period, but such a long modulation period has not been 
found yet (see Sec l3.2t . Analysis of the residuals also did not show 
any signs of modulation. 

Only one single RV measurement for RZ Cet is available 


RZ Cet model 2 
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RZ Cet model 2 - parabola subtracted 



Figure 5. O—C diagram of RZ Cet together with model 2 (top panel) and 
variation after subtracting the parabolic trend (bottom panel). Symbols are 
the same as in Fig.[T] 


dColacevichll 195(]l) . which is inappropriate for our purposes. Spec¬ 
troscopic measurements are definitely highly needed to confirm the 
binarity of RZ Cet. 


5.4 S Com 

S Com is another promising candidate for a binary system. The 
information about possible cyclic changes with a long per iod was 
mentioned without closer details in lLe Borgne et al.] OOOTh . Proba¬ 
ble LiTE with a period of about 106yr is well demonstrated in our 
O—C diagram (model 1, the top panel in the Eig.[^. Currently, the 
pulsation period is at its the shortest, but in the next few years it 
should start rising again. 

Our second model of LiTE of S Com including the parabolic 
trend (the middle panel in Eig. |6j gives a shorter orbital period of 
90 yr and also lower orbital parameters such as the projection of 
the semi-major axis or the limiting mass of the secondary compo¬ 
nent. The rate of period changes = —l.tCW x 10^** dd^* 
significantly differs fr om the value —26.4(1.1) x 10^** dd^* 
jLe Borgne et al.l [2007h mainly due to inclusion of LiTE in the 
model. 

There are two RV datasets available in the literature (Table|4]l. 
However, expected changes of the systemic velocity caused by 
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Figure 6. O—C diagram of S Com together with our model 1 (top panel), 
model 2 (middle panel) and variation after subtracting the parabolic trend 
(bottom panel). Symbols are the same as in Fig.[T] 


binarity with a semi-amplitude of only 1 . 6 kms^* (model 1 ) or 
1.3 km (model 2) are completely saturated by the RV changes 
caused by pulsations. 


Figure 7. O—C diagram of SU Dra together with our model 2 (top panel) 
and variation after subtracting the parabolic trend (bottom panel). Symbols 
are the same as in Fig.[T] 


5.5 SU Dra 

0—C variations of the star SU Dra cannot be explained by a sim¬ 
ple parabolic trend. As in other mentioned cases, we proposed that 
period changes are caused by the presence of an unseen body. Ac¬ 
cording to our model, combining parabola and LITE (Fig. |3, the 
low-mass companion (9Jt2,min ~ 0.09 DTIq) should orbit SU Dra 
with a period of 70.7 years. The parabolic shape of 0—C sug¬ 
gests a period change rate of 16 . 6 )')j '2 x 10 ^** dd^* which is al- 
most the same value as 17.0(1.0)x 10^*^ dd^* dUe Borgne et al.l 
l2007lf'^ An additional scenario (simple LITE) was tested and sim¬ 
ilar to other proposed systems, the found orbital period of about 
149 yr is very uncertain since only an incomplete cycle is available 
in the 0—C diagram. 

SU Dra, as the brightest star from our sample, has the largest 
number of sources with available RV measurements (Table|4j. Un¬ 
fortunately several values are dubious and the remaining data seem 
to contradict a binary hypothesis. However, the semi-amplitude of 
predicted RV changes Ki is quite small - only 0.85 km s^* which 
is hardly detectable. Therefore this test should not be considered as 
reliable. 


He calculated also third-order fit for SU Dra with 13.4x 10 **dd *. 
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Figure 8. Light curve of VX Her in V band obtained in order to detect 
eclipses at MUO during 13 nights (April - August 2014) together with 
the polynomial model (top). The residual after model subtraction (bottom) 
shows no evident variations larger than 0.02 mag. 


5.6 VX Her 

VX Her was proposed t o be an eclipsing binary by 
iFitch. Wisniewski & Johns^ i ll 9661) - see the discussion in 
Sect. 12.11 for more details. Unfortunately, no other observation 
from literature nor our observations (Fig.[^ confirm Fitch’s results. 

The shape of the 0—C diagram (the top panel of the Fig.|9) of 
VX Her is apparently a parabola with a small systematic deforma¬ 
tion. The residuum after subtraction of the parabolic trend (bottom 
panel of the Fig. |9} contains a small wave which can be described 
by a model of LITE with a period of about 86 yr. Nevertheless, 
the existence of LiTE is a bit questionable due to the high scatter in 
older maxima (photographic, visual). At least our result can be con¬ 
sidered as a limiting value for LiTE in VX Her (the limiting mass 
of a second component could only be 0.05 DJIq). The value found 
for the period change rate —42.71^®® x 10^** dd^* is similar to 
the value —40.5(9) x 10^** dd^* from iLe Borgne et ^ ( |2007|) . 

All the 0—C variations can also be described using only sim¬ 
ple LiTE with a longer period. In such a highly unlikely scenario 
the possible orbital period would be about 175 yr, and the limiting 
mass of the secon dary star would th en b e about 0 . 8 pyip) . 

According to lAzarnw i d 19631) and lWund^ d 19901) the pulsa¬ 
tions are modulated with a Blazko period of 455.37 d with quite a 
high amplitude of 0—C changes of 0.013 d. This periodicity was 
found after subtraction of the parabola in the 0—C diagram with¬ 
out verification of light variations. Later results from Wunder are 
based only on photographic, visual and 4 photoelectric measure¬ 
ments of maxima timings and they should be verified. Moreover 
these changes may not be caused by the Blazko effect, but by 
LiTE. In our residuals after subtracting the parabola, which are 
based on photoelectric and CCD measurements, no indications of 
modulation are apparent. Since these residuals have a standard de¬ 
viation of 0.0015 d, changes in amplitude of a magnitude larger 
should definitely be detectable. Erequency analysis of residuals 
from ph otographic and visual measurements (including measure¬ 
ments of lWundeilll99(lh show no indications of periodicity . In ad¬ 
dition , Blazko modulation was undetectable in ASAS data dSkar 
I 2 OI 4 I) and in our measurements (Eig.[8]l. Therefore it is very likely 
that VX Her is currently unmodulated. 

RVs from two available sources (see Table 0 have insuffi¬ 
cient quality to confirm the binarity, because the range of RV vari- 
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Figure 9. O—C diagram of VX Her together with our model 2 (top panel) 
and variation after subtracting the parabolic trend (bottom panel). Symbols 
are the same as in Fig.[T] 



VX Her model 2 



ations is from —410 to —340 km s *, while our model suggesting a 
O.O 5 -OJI 0 companion predicts the semi-amplitude of RV to be only 
0.8 km s-'. 


5.7 SS Leo 

The O—C diagram of SS Leo shows large changes in the pulsa¬ 
tion period (Fig. that can be explained by long-period LiTE. 
Apparently it could not be explained as secular period change with 
additional low-amplitude LiTE. We found a probable orbital period 
with length of about 111 years. Therefore, only one cycle was com¬ 
pleted since the time of the first observation. SS Leo is known as a 
regular pulsating star without the Blazko effect dSkar OI 2 OI 3 I) . and 
we did not either notice signs of the Blazko effect. The proposed 
companion should have a minimum mass of 0.38 9Jl©. 

Measurements of RV were found in four publications (see Ta¬ 
ble 0. Unfortunately, three of them contain values obtained over 
4 years, thus practically in the same phase of the supposed orbit. 
During this time the difference in systemic velocity based on our 
LiTE model should change only within 0.16 kms^*, and this dif¬ 
ference is below the accuracy of these RV measurements. Never¬ 
theless, the pulsation phase is well covered by RVs and values vary 
in the rang e of 120-205 kms^*. One of the values published by 
Ia^ d 19731) (-20kms-') is very different from the others and is 
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Figure 10. O—C diagram of SS Leo together with our model 1. Symbols 
are the same as in Fig.[T] 


probably incorrect. Again it is impossible to confirm the binarity of 
SS Leo using available RV measurements. 

However, in May 2015 the systemic velocity should be lower 
by about 4kms^* than for RV measurements from the last three 
sources. We propose SS Leo as an easy target for confirmation of 
the binary nature by RV measurements. 


5.8 AV Peg 

AV Peg’s 0—C diagram has a parabolic shape (Fig. [IB. In 
the residual, after applying a qua dratic ephemeris, a wave was 
detected, which was proposed by ISzeidl. Olah & MizseJ l ll986ll 
as a sign of LiTE. Their preliminary value of the orbital pe¬ 
riod was estimated to be 64 years. Our approach, in which pa¬ 
rameters of a quadratic ephemeris and LiTE are fitted simul¬ 
taneously (top panel of Fig. dB, yields a shorter orbital pe¬ 
riod of 48 yr and the rate of secular change of the pulsation pe¬ 
riod Fpuis = 4-45.96^4, X 10^" dd^* which is very similar to 
values -f47.3(8) x 10~‘^ dd~^ jSzeidl. Olah & MizseJfl^ and 
4-45.8(5) X 10^** dd^* jLe Borgne et al.ll2007h . The ^fference in 
the length of the orbital period can be explained by using about a 
30 year long er dataset with more accurate C CD measurements than 
was used bv ISzeidl. Olah & Mizsej ( Il986h . Our model of LiTE in 
AV Peg gives the second lowest minimum mass of the secondary 
component of all stars in our sample (only 0.06 StU©). When apply¬ 
ing the model only to CCD and photoelectric measurements and 
subtracting the parabolic trend (bottom panel of the Fig. dB, the 
LiTE is more pronounced than in the case of all measurements 
(middle panel of the Eig. dB. 

RVs found in 3 sources (Table II were obtained in 
phases where the systemic velocity should differ by only about 
0.26 km s^*. Thus the variations in RV due to binary motion are 
undetectable in these measurements. The range of RVs due to pul¬ 
sation of the staiF^ is from —93 to —31 kms^*. In any case, the 
semi-amplitude of RVs caused by the proposed orbit is too low 
(only 0.7 km s^') and will be very difficult to confirm. 
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One value from lAbl|jl973h substantially differs from the others. 
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Figure 11. 0—C diagram of AV Peg together with our model 2 (top panel). 
The middle panel demonstrating LiTE shows all data after subtracting the 
parabola, the bottom panel shows the same only for CCD and photoelectric 
measurements. Symbols are the same as in Fig.[T] 


5.9 AT Ser 

More than one orbital cycle is covered for times of maximum for 
AT Ser (Fig. ll2L From the amplitude of changes in 0—C it is evi¬ 
dent that the possible companion has a high mass. Our model, with 
an orbital period of 86 years, predicts a minimum mass of the sec¬ 
ond component of 1.9 QJl©. Therefore the companion should be a 
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Figure 12. The O—C diagram of AT Ser together with our model 1. Sym¬ 
bols are the same as in Fig.[T] 

degenerate remnant in the form of a neutron star, or more likely in 
the form of a black hole. 

RVs were found in 3 sources (Table|4ll. Unfortunately, these 
9 values do not even cover a complete pulsation phase. Measured 
RVs vary between —100 and —40 km s^*, but the full range of vari¬ 
ations in velocity will be a bit larger. The calculated semi-amplitude 
of the orbital RV curve for this system is one of the largest - about 
7.8 kms^^. Currently (May 2015), the systematic offset in RVs 
should have t he highes t differ ence of about 15.5 km s^* other than 
for RVs from Ipemlevl ( Il993h . AT Ser is, therefore, another easy 
target for proving binarity using RVs. 
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5.10 RV UMa 

RV UMa is a well known modulated star with secular period 
changes accompanied with LiTE covering more than one proposed 
orbital cycle (Fig. [nil. When the 0—C dependence is considered 
as LiTE with a period of 66.9 yr supplemented by a slow evo¬ 
lutionary period change (parabolic trend with the rate of Rpyis = 
4.7 X 10^** dd^*, the companion should have a minimum mass of 
0.25 OTq. 

Except for these large-amplitude period changes, RV UMa 
shows additional sm all changes cau sed by the Blazko effect with 
interesting behaviour. lKanvdr ( ll976h analysed the period behaviour 
of RV UMa and found that the Blazko period (about 90 d) changes 
in the antiphase with the pulsation period. A more detailed analy¬ 
sis of peri od evolution and th e Blazko effect, performed bv iHurtal 
( l2007h and iHurta et al.l ( l2008l) . confirmed Kanyo’s results and com¬ 
plemented it with finding that between 1946 and 1975 the changes 
in th e pulsation and m odulation periods were parallel. In addi¬ 
tion, iHurta et ^ (l2008h gives the range of the modulation-period 
changes between 89.9 and 90.6 dB Nevertheless, none of these 
authors discussed the possible binary nature of the cyclic varia¬ 
tions of the 0—C diagram shown in the Fig. [T^ After removing 
LiTE and the parabolic trend we searched for periodicity caused 
by the Blazko effect. We found a modulation period of 89.9(1) d 
(Fig.Ql. The high scatter probably corresponds to the variation in 
the length of the Blazko period. 

iKanvdl il97d) found that in 1961 the Blazko period was 92.12 d. 


Figure 13. O—C diagram of RV UMa together with our model 2 (top panel) 
and variation after subtracting the parabolic trend (bottom panel). Symbols 
are the same as in Fig.[T] 


We attempted to confirm variations in the O—C diagram by 
RV variations. Three datasets of RV measurements that cover the 
pulsation cycle well were found in literature (see Table 0. The 
measured values of RV were compared with RV values corrected 
for systemic velocity caused by orbital motion (parameters from 
our model of LiTE). Unfortunately, confirmation failed due to a 
small semi-amplitude of K\ of 2 .2kms~* and due to ch anges in 
RV during the Blazko effect - IPreston & Sninr^ ( Il967h clearly 
presented evolution of the RV curve for RV UMa in different 
Blazko phases (change of the shape and amplitude of more than 
20 kms-'). 

_ Our r esults, together with those from iKan 3 ( Il976l) . iHurtal 

( l2007h and ljurcsik et all ( l2002h . suggest that there could be some 
connection between orbit al motion and variations in pulsation and 
Blazko periods. Although iHurtal ( l2007h did not deal with binarity, 
67-yr antiparallel pulsation and modulation period changes, appar¬ 
ent from his Fig. 4, indicate this possibility. Similar but parallel 
behaviour was observed in XZ Dra jjurcsik et al.l2002h . and we as¬ 
sume that the connection with bin arity should also be in vestigated 
in other Blazko stars, e.g. XZ Cyg jLaiciuvze et al.ll2004l) . RW Dra 
(lFirmanvuklll978h . 
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Figure 14. 0—C residuals of RV UMa phased with the modulation period 
of 89.9(1) d. 

5.11 BB Vir 

iKinman & Carretti ( Il992h supposed BB Vir to be a member of a 
binary with a hot HB star (already discussed in Sect. I2.5l l. When 
we check Fig. [15] we immediately see that the 0—C dependence is 
deformed and could not be fitted with a simple parabola. One of the 
possible explanations is that the 0—C is a part of the LiTE with an 
incomplete cycle. Assuming this less probable scenario, we found 
that the companion with a minimum mass of about 0.95 DJI© should 
orbit around the RR Lyrae component at a minimum distance of 
21.6 au with an orbital period of about 168 years. The relatively 
high minimum mass suggests a degenerate remnant. 

Nevertheless, the second model combining the parabolic trend 
with LiTE (Eig. fTst is supposed to be more probable. We found 
a period of about 93 yr and the rate of period changes Rpui, = 
+ 36 . 3 ('(q 3 X 10^^* dd^*. The obtained period-change rate dif¬ 
fers from value +46.0(3.8) x 10^*^ dd^* ^Le Borgne et alj[2007h 
which can be explained by adding LiTE to our model. 

Because BB Vir is known as a star with modulation in the 
GCVS catalogue, we unsuccessfully searched for additional pe¬ 
riodicity in residuals. In addition, the modulation period is not 
known, and n o signs of modulation were found in ASAS data by 
[Sk 3 j2014t) . Therefore the Blazko effect in BB Vir is at least 

questionable. _ 

Only one single RV measurement —15.6 km s" l lAblll973h is 
available in literature, but it is insufficient to tell us anything about 
the binarity of BB Vir. 



5.12 Other RR Lyrae stars with long-term cyclic variations 
in 0—C mimicking LiTE 

There could be several other RR Lyrae type stars that exhibit cyclic 
changes that could possibly be interpreted through LiTE. Never¬ 
theless, the shape of variations often changes from cycle to cy¬ 
cle which is difficult to explain by simple LiTE. A typical exam¬ 
ple is the star RR Lyr itself. Its semi-regular variations after the 
gap at about JD 2435000 resemble the LiTE with a 14.1-yr period 
(Fig.QDi and a companion with a minimum mass of 2 DJI© with 
an eccentric orbit with e ~ 0.4. Unfortunately, the general trend in 
period evolution of RR Lyrae is more complex with the gap, and 
different slopes before and after the gap. The situation is similar 
with AQ Lyr (Fig. [iTt which underwent two cycles during the in- 
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Figure 15. O—C diagram of BB Vir together with our model 2 (top panel) 
and variation after subtracting the parabolic trend (bottom panel). Symbols 
are the same as in Fig.[T] 


terval JD 2424000 -2438000 with a period of 19.8 yr well approx¬ 
imated by a highly eccentric orbit (e 0.75 — 0.95) of a massive 
component 6-9 DJI©. Before and after this interval the star behaved 
differently. 

Another candidate for binarity is AE Peg showing high ampli¬ 
tude variations in O—C (Eig.ll8t with a period of about 99 years, 
e ~ 0.3, and a semi-major axis ai sini ~ 56 au. The secondary 
component should be a massive black hole with a mass of about 
18.6 DJI©. Maybe a binary explanation is not correct for these ob¬ 
jects, but hypotheses employing other effects discussed in Sect. 
are equally uncertain. 

These three mentioned objects have known Blazko modula¬ 
tions with lengths of tens of days. Detected long-term cyclic varia¬ 
tions can be a manifestation of multiple Blazko effects. These ob¬ 
jects are good examples of stars with complex behaviour that looks 
like LiTE in some parts, but in others are undetectable. In the case 
of AQ Lyr and RR Lyr, a time-base longer than 2 and 4 cycles, 
respectively, was necessary to detect a completely different kind 
of variations in the O—C diagram. Additional observations will be 
highly appreciated. 
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Figure 16. O—C diagram of RR Lyrae (after the gap at about JD 2435000) 
shows also cyclic changes which could be described by the LiTE model. 
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Figure 17. O—C diagram of AQ Lyrae that shows cyclic changes during 
a short interval (JD 2424000-2438000) which could be described by the 
LiTE model. 
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6 DISCUSSION 

Cyclic changes in 0—C diagrams can be interpreted in various 
ways. The time scale of the changes that we studied is crucial. 
Thus far the longest observed process in RR Lyrae stars causing 
cyclic changes (not necessarily strictly periodic) in the 0—C di¬ 
agram, is the Blazko effect. Observed periods are maximally in 
the order of years, while we analysed period-evolution in the order 
of several decades. Another possibility to explain cyclic long-term 
period variations could be the pr esence of a sola r -like magnetic 
cycle. This idea was proposed by IStothersI dlQSOl) . IPerekas et al.l 
( 120041) argued that the relative change in radius during the period- 
change cycle of BE Dor is in agreement with the rel ative change 
of solar radius during its magnetic cycle dNoe l l2004h . and there¬ 
fore a hydromagnetic interpretation could be correct. The problem 
is th at a dipole magnetic field in RR Lyrae stars was not detected 
(e.g.. lChadid et al.ir2004lKolenberg & Bagnuldr2009t) . and it is not 
clear whether solar-like magnetism is even possible in pulsating 
HB stars. Perhaps some unknown physical mechanism could be re¬ 
sponsible for the long-term cyclic evolution of the pulsation period. 
However, the most simple, and currently most probable explanation 
is through binarity, which we assume in this paper. 

Our target selection gives preference to orbits with wide tra¬ 
jectories and with long orbital periods lasting for many decades. 
Such systems are, therefore, definitely well detached without mass 
transfer between components taking place during their life cycles. 
This should allow direct determination of the true RR Lyrae mass 
which is not influenced by stellar evolution in a close binary system 
- a chance keenly awaited by all theoreticians dealing with stellar 
evolution and pulsations. However, wide trajectories make eclipses 
highly improbable, not even mentioning eclipse prediction, which 
strongly degrades the chance of precisely determining the mass of 
an RR Lyrae component in this way. Observation of eclipses is 
therefore a very un likely event. Maybe the dimming o f the light 
of VX Her noted by lpitch. Wisniewski & JohnsonI (Il966h was such 
a case. In spite of the slim chance for eclipse detection, wide tra¬ 
jectories can be used with the advantage of determining RR Lyrae 
mass on the basis of visual binary motion together with RV mea¬ 
surements. _ _ 

Howev er, _ Ae recent studies by iLi & Oia^ (12014l) and 

iHaidu et"^ ( 12015l) detected systems with much shorter orbits only 
in the order of yearj^ giving a better chance for catching an 
eclipse. OGLE and Kepler measurements span maximally a few 
years. Thus long-period binaries remained hidden in these observa¬ 
tions. A great advantage of the data from OGLE and space surveys 
against our sample is that they observed many thousands of targets 
almost continuously, while we analysed bright field targets with a 
large amplitude of variations in their 0—C, that are sparsely ob¬ 
served separately by different observers, and not regularly and with 
similar cade nce as in the case of K epler or OGLE. The time-span of 
the data that iLi & OianI ( 12014l) and iHaidu et al.l ( 12015h used, caused 
that their results are strongly biased towards short orbital periods. 

Nevertheless, a project established for period studying of RR 
Lyraes, using their maxima timings regularly and densely observed 
by the TAROT telescopes dLe Borgne et aljEoo4 l2007h . can be 
considered as comparable to the OGLE. Detailed analysis of all 
0—C diagrams (not only a visual selection as was done in this pa- 


Figure 18. 0—C diagram of AE Peg that shows probable cyclic changes 
which could be described by the LiTE model. 


Binary system TU UMa iLiska et ^|2Q15|) with a 23.4yr period is in 
the middle of the orbital period interval of the mentioned survey systems 
and those from this paper. 
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per) for these TAROT objects would reveal systems with periods of 
several years. 

Another problem which relates to the length of the proposed 
orbital periods, are the uncertainties surrounding the choice for 
modelling of the O —Cs (with or without secular parabolic trends). 
iHaidu et^ l l2015h applied secular changes together with LiTE by 
default. Generally, our models with parabolic trend give much 
shorter periods than models with LiTE alone. 

Companions of all stars in our sample should be objects with 
minimum mass larger than about 0.05 OT©. In the case of VX Her 
and AV Peg, their companions could be brown dwarfs when the 
orbits would have an inclination close to 90 degrees. In other ex¬ 
amples the companions are probably low-luminous low-mass main 
sequence red dwarfs, or stellar remnants. This corresponds well 
with the absence of information about the peculiar colour of any 
of the stars in literaturjn. If the companion was a luminous star of 
a different colour than the RR Lyrae component, we would observe 
significant colour excess in comparison with the average colour of 
RR Lyraes. 

In RZ Cet and AT Ser our models suggest high-mass compan¬ 
ions in the form of a degenerate remnant with a minimum mass 
of about 1.15 9)10 and 1.90 ORq, respectively. Progenitors of these 
massive remnants produced heavier elements without a doubt and 
definitely should enrich their RR Lyrae partners with metals - as it 
was proposed bv iKennedv et al.l l l2014h . Thus observation of higher 
metal abundances should be an independent and relatively easy way 
to confirm the binary nature of RZ Cet and AT Ser. 

In addition, the wide orbits of our systems represent a perfect 
opportunity to prove their binarity through direct imaging. The dis¬ 
tance, calculated using the observed and average absolute magni¬ 
tude (0.6 mag|3 together with a computed semi-major axis, gives 
an angular distance between binary components ranging from 1 to 
13 mas for our sample stars. This is ab out ten times more than for 
stars in the Galactic bulge studied by iHaidu et"^ l l2015h . In ad¬ 
dition, the appropriate orientation of the orbit can further double 
these numbers. On the other hand, the angular distance could be 
much lower. The resolution of a few mas should be reachable by 
state-of-the-art ground-based and space instruments. Because the 
proposed companions are low-mass, low-luminous red stars, the 
detection would be less demanding in the red part of the spectrum 
where the difference between members of the system is less. How¬ 
ever, even in IR the difference would be several magnitudes which 
could cause some problems in detection. The most promising can¬ 
didates with the largest semi-major axes are AT Ser (13 mas), RS 
Boo (9 mas, model 1) and RZ Cet (8 mas). 

We did not notice any correlation between the presence of the 
Blazko effect and binarity. However, some unknown mechanism 
connected with binarity can induce variation in the length of the 
Blazko period as was discussed for RV UMa in Sect. 15.101 Also 
metallicity probably does not influence the fraction of RR Lyraes 
in binary stars. 


7 SUMMARY 

In this study we focused on RR Lyrae type stars with suspicious 
0—C diagrams showing cyclic variations or some indications of 
possible long-term periodic changes, preferably of large amplitude. 

iKinman & Carrett j ll992h found BB Vir to be somewhat bluer than 
other RR Lyrae stars. 

** Hipparcos parallaxes are unreliable for the studied stars. 


We attribute them to the consequence of binarity manifesting itself 
as the Light Time (Travel) Effect. We attempt to provide an exten¬ 
sive overview of known, or more accurately, suspected RR Lyrae 
type pulsators in binaries, and we discuss problems involving ob¬ 
servational and analytical methods for the confirmation of the bi¬ 
narity. We also widely discuss various effects which can influence 
the appearance of 0—C diagrams to establish a solid basis for in¬ 
terpretation and strengthening of our results. 

We proceeded from the collect ion of maxima timings provided 
in the GEOS RR Lyrae database dLe Borgne et al.ll2007ll and se¬ 
lected several stars with 0—C diagrams whose appearance could 
be interpreted as the consequence of binarity. In addition to data 
from the GEOS database, we determined maxima timings for se¬ 
lected targets from large sky surveys and our observations. Eor the 
modelling of the LiTE, we used recently developed co de based on 
the no n-linear LSM and bootstrap-resampling method dLiska et al.l 
l2015h . LiTE for eleven stars were modelled and possible orbital pa¬ 
rameters were calculated (Table 0. Orbital periods of hypothetical 
binaries range from 47 (AV Peg) to 147 years in RS Boo (model 1). 
The projection of the semi-major axis of the pulsating component 
a I sini was estimated as 1 au in AV Peg and about 20 au in AT Ser. 
Minimum masses of companions were obtained according to our 
expectation (see Sect.[Tll. They were either very small (hundredths 
to tenths of D)!©) or higher than the solar mass. All proposed sys¬ 
tems have eccentric orbits, but none of them has e > 0.9. The ma¬ 
jority of studied stars has an eccentricity between 0.1 and 0.5. Only 
RS Boo and VX Her could have more eccentric orbits. 

Although in many targets we could not clearly decide which 
model of 0—C variations is correct (only LiTE, or LiTE super¬ 
imposed on a parabola), our models provide a prediction for the 
period evolution in future years, as well as for RV which could 
be expected to be observed in years to come. Consequently only 
the future can provide a reliable answer to the question whether 
RR Lyrae stars analysed in this study are really bound in binary 
systems or not. Nevertheless, our study substantially extends the 
group of suspected RR Lyrae variables in binary systems with pre¬ 
liminary, estimated orbital parameters. 
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